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On Abstract grad—div Systems. 
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Abstract. For a large class of dynamical problems from mathematical physics 

/ 0 —C* \ 

the skew-selfadjointness of a spatial operator of the form A = I ^ ^ j, where 

C : D {C) C Hq —)■ Hi is a closed densely defined linear operator, is a typical 

/ 

property. Guided by the standard example, where C = grad = 


(and 


\dn 


—C* = div, subject to suitable boundary constraints), an abstract class of op- 

/c. , 

is introduced (hence the title). As a particular application 


erators C = 


\Cn 


we consider a non-standard coupling mechanism and the incorporation of diffusive 
boundary conditions both modeled by setting associated with a skew-self adjoint 
spatial operator A. 


Keywords and phrases: Evolutionary equations, Gelfand triples, Guyer-Krumhansl heat conduction, Dy¬ 
namic boundary conditions, Leontovich boundary condition 
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0 Introduction 


In a number of studies |211 EHl EH 113 EHl |22] it has been illustrated, that typical initial 
boundary value problems of mathematical physics can be represented in the general form 


(doM +A)U = F, 

where A is skew-selfadjoint, indeed commonly of the specific block matrix form 


A = 


0 -C* 
C 0 


( 0 . 1 ) 


( 0 . 2 ) 


with C : Xi C Xq Y a. closed densely defined linear operator between Hilbert spaces Xq 
and Y with Xi = D(C), see e.g. |15[ 116) . The operator Ai is referred to as the material law 
operator, which in the situation discussed here is a linear operator acting on a Hilbert space 
realizing the space-time the problems are formulated in, m isi EH ES]. 

The main purpose of this paper is to focus on the operator C in this construction of the 
operator A : D (C) x D (C*) C Xo©T —)• XqTT when Y is itself a direct sum of Hilbert spaces. 
In such a situation we shall loosely refer to a system of the form m as an abstract grad-div 
system. The guiding example, which at the same time motivates the name, is to take for C : 

Xi C Xq —>■ Y the differential operator V = I <92 I with domain Xi = Hi (H), the Sobolev 

space of weakly differentiable functions with L^-derivatives and vanishing boundary data. 
The spaces Xq and Y would be in this case (H) and (= (H) © (H) © (H)), 

respectively. The corresponding equation of the form (10.11) would lead in particular to a model 
for the acoustic wave propagation or - depending on the material law operator Ai - for the 
dissipation of heat with so-called Dirichlet boundary condition, see e.g. ns [231 |32]. The 
adjoint operator is the negative weak divergence C* = — div. It is 


C 



(0.3) 


The idea in this paper is to replace the role of the partial derivatives in ()0.3p by general opera¬ 
tors in general Hilbert spaces, hence the term abstract grad-div systems for the corresponding 
evolutionary systems associated with the skew-selfadjoint operator A constructed according to 
dOj). That the study of this particular class of skew-selfadjoint operators leads to interesting 
applications is illustrated by three examples. First we consider the so-called Guyer-Krumhansl 
model of thermo-dynamics, ElEllH], as a particular instance of this construction. The last 
two examples illustrate that the concept of grad-div systems is well-suited to study dynamic 
boundary conditions. A first implementation, leading up to the concept of grad-div systems 
introduced here, has been established in the context of a particular class of boundary con¬ 
trol problems, |18[ EZ]. Our investigation of grad-div systems has been motivated by the 
discussion of a heat conduction problem with a dissipative dynamic boundary condition in 
@]. The generalization to the context of evolutionary equations is the topic of our second 
example. The last application example is connected to the Leontovich boundary condition 
of electrodynamics, see e.g. |121 (S) [25]. We shall discuss two different implementations. The 
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first one is based on classical bonndary trace concepts (see e.g. mm). the second approach 
nses a bonndary data concept developed in IE1[IT| with extensive nse in |28) . which has the 
advantage that no constraints on the qnality of the bonndary are incnrred. 

In Section [1] we start by presenting the fnndamental constrnction of abstract grad-div sys¬ 
tems. The remaining section is then dedicated to illnstrating the nsefnlness of the concept by 
examples. 


1 Construction of Abstract grad—div Systems. 


In this section, we shall reconsider the concept of the adjoint operator of a densely defined, 
closed linear operator C, specifically in order to deal with the fact that the image space Y of 
the operator C is given as an orthogonal snm of Hilbert spaces. Let ns first provide a precise 
definition of what we wonld like to call an abstract grad-div-system. 

Definition 1.1. Let C : Xi C Xq —>■ T, be a densely defined, closed linear operator with 
domain Xi between Hilbert spaces Xo,Y. We shall refer to a system of the form (10.11) with 
A generated via dEl, as an abstract grad-div system, if Y given as a direct snm, i.e. Y := 
0fcG{i n} for Hilbert spaces Y^, A: G {1,... , n}, n G N. 


As a matter of jargon we shall say that the abstract grad-div system is generated by C. If tYk 
denotes the canonical isometric embedding of Yfc into Y then, with Ck ■= /c G {1,... ,n}, 

we have 


Cx = Cix © • • • © CnX = 



( c A 


\cj 


X G 



= Y 


for X G Xi. 


The aim of this section is to give a characterization of the adjoint C* of C. However, before 
stating and proving the respective theorem, we shall explore as a first elementary example the 
name-giving case mentioned in the introdnction: 


Example 1.2. Let n G N, H C R” non-empty and open, Xq := (H) and Xi := iLi(H), 

the Sobolev space of weakly differentiable fnnctions in L^(H) with vanishing Dirichlet trace, 
i.e., the closnre of compactly snpported continnonsly differentiable fnnctions with respect 
to the Fi-norm / ^ Then clearly, with Y = (H) 

the operator C '■= grad : Hi{0) C (H) —?■ (H), defined as the closnre of the 

classical vector-analytical operation grad restricted to smooth fnnctions with compact snpport 
in H, generates an abstract grad-div-system. In the same way, choosing Xi = Lfi(H) instead, 
we get that C = grad also generates an abstract grad-div-system. We recall that the difference 
between the operators grad and grad lies in the additionally prescribed homogeneons Dirichlet 
bonndary condition for elements in the domain of grad. 


To clarify notation, we need the following definition. 
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Definition 1.3. Let X, V be Hilbert spaces and L : X ^ Y a bounded linear operator. Then 
we denote the dual operator of L by L' ; T' —)• X' defined by 

{L'y') {x) := y'{Lx) {y' eY',x e X). 

Moreover, we may consider a densely defined linear operator S : Xi C Xq —>■ Y such that 

Ls-.Xi^Y 
X ^ Sx 

is a continuous linear operator {S need not be closable) and define the operator : Y ^ X[ 
by := L'g o Ry. Here Ry' : Y ^ Y' denotes the Riesz isomorphism, given by 

{Ryy) (z) := (ylz)r (y, z G Y). 

Note that we have by definition 

(S'^y) (x) = (LsiRv'y)) (x) = (Ryy) (Sx) = {y\Sx)Y 
for y G y, X G X. As a matter of convenience this can be written suggestively as 

{S^y\x)x^ = {y\Sx)Y, 

where ( • | ■ )xq denotes here the continuous extension of the inner product of Xq to a duality 

Rx^ 

pairing on X[ x Xi. Note that Xi ^ Xq =° Aq X[ is a Gelfand triple. 

Lemma 1.4. Let Xo,Xi,Y be Hilbert spaces, Xi C Xq dense, C: Xi C Xq Y a closed 
linear operator. Then the adjoint C* of the operator C is densely defined and its adjoint C* 
is given by 

c* = c^r\{Y® Xq) , 

which is the same as to say 

C* = {{y,x) £ Y ® X[ I C^y = X and x G Xq } . 

Proof. Note that, by definition, C* C and C* ®Y 0 Aq. Hence, C* C (Y © Aq). For 
the remaining implication, let (y, x) G A 0 Aq. Then we have 

(y, x) G C* V(/> G Ai: {y\C4>)Y = {x\fi)xo 
V(^ G Ai: C'^y((/)) = (x|(/))xo- 

This establishes the assertion. □ 


The latter observation has been employed in more concrete situation in control theory, see 
e.g. |171 ITHl [19] and, in particular, the references in )19) . 


Lemma 1.5. Let C : Ai C An —)• Y generate an abstract grad-div system with Y = 
Then 

/ yi\ 


C^={Ct--- C:):Y3 


I-)- 




\yn ) 


k=l 


7 




Proof. We have for x £ Xi, y = 


/ yi\ 


\ynj 


= 2/1 


yn^y, 


(C'^y)(x) = {y\Cx)y = ^ {yk\Ckx)Y^ = {Ctyk) {x ), 

k=l k=\ 

hence, = Y2=i C^yk- 

The latter two lemmas yield the proof of the following theorem: 


□ 


Theorem 1.6. Let C generate an abstraet grad-div system with C = 


Cl 

Cn 


Then 


C* = {Cl ■■■ Cl)r^{Y®Xo) 

= \ {{yi,---,yn),x) G y©Xo 


X = 


Y^tyk^Xo 


k=l 


Let us apply the latter theorem to the case C = grad, as it was defined in Example 11.21 

Example 1.7. Let Xq = L‘^{Ll),Xi := Hi{Ll) and Yk '■= Lp‘{Pl) for k G {1, ... ,n} and some 
non-empty, open set Q, C Recall that C = grad generates an abstract grad-div-system. 

Cl ^ 

with Ck '■ Hi{kl) C (n) — Lp‘{Q) is defined as Ckf = dkf for 


As in this case, C = 


C„ 


f G Hi{Ll), we obtain in particular that 

(Ctg) if) = (t/|4/)L2(o) = - {dkg) (/) 

for g G L^(fl),/ G Coo(fl), where here dkg is meant in the sense of distributions. Here 
(7oo(f^) denotes the set of arbitrarily differentiable functions with compact support in In 
consequence 

/ gi\ 

~ ^ ^ dkgk 


C" 


\gnj 


k=l 


and thus, C* = — div is the negative distributional divergence on L?‘{Ll) vector-helds restricted 
to those, whose divergence is representable as an L^(H)-function. 

An immediate application of Theorem 11.61 is the following corollary, which will be of interest 
in the next section: 

Corollary 1.8. LetC generate an abstract gxa.d-di'v-system. Assume that there exists a closed, 
densely defined operator Ci with 

/Ci\ 

0 


C 


V 0 / 


Cl 

Cn 


= C. 
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Then 


C* = (Ci* 0 • • • o) 

Proof. It suffices to observe that C* C 0 • • • 0 


C 0 • • • o). 

) and to apply Theorem 11.61 


□ 


The next statement contains a typical situation of operators of the form C 
rise to the generation of abstract grad-div systems. 




CnJ 


giving 


Proposition 1.9. Let Xq, Yq, ... ,Yn be Hilbert spaces, Cq : D{Cq) C Xq —)■ Tq densely defined 

and closed. Denote Xi := ^Il(C'o), \/| • p + \Cq ■ and let Ck G L(Xi,Yk), k G 

Then 


C 


/Co\ 

: : DiCo) C Xq ^ 0 ^ 

\Cn/ fce{0,...,n} 


generates an abstract grad-div system. 


Proof. As Co is densely defined, so is C. Thus, closedness of C is the only thing remaining to 
check. For this let be a sequence in D(Cq) convergent in Xq its limit being denoted 

by X G Xq and with the property that ■= limj_>.cxD CkXj exists in Yj,, j G {0,..., n}. By the 
closedness of Cq, we infer x G D[Cq) and Cqx = yQ. Thus, Xj ^ x in Ai as j —)• oo. The 
continuity of for all k G {1,..., n} yields the assertion. □ 

Remark 1.10. Note that the operators Ck, k G {1,... ,n}, need not be closable with domain 
D{Cq). The following example illustrates this fact. Take Xq := L^(0,1), Cq := d: H^[0, 1) C 
L^(0,1) —)• T^(0,1), / e-)• f with f denoting the distributional derivative. Then Cq is densely 
defined and closed. It is known that Ffi(0,1) C C[0,1] (Sobolev embedding theorem). Hence, 
Cl/ := (5 {i/ 2}/ •= /(1/2) for / G Ffi(0,1) defines a continuous linear operator from Ffi(0,1) 
to C. It is easy to see that Ci: H^i(0,1) C L^(0,1) —)• T^(0,1) is not closable. 


2 Some Applications 


Before we come to specific applications, we introduce some classical differential operators, 
which we will need in their description. Throughout, let H be a non-empty open subset of M"', 
n G N. 

Definition 2.1. We define the operator grad as the closure of the operator 

Coo(H) C L^{n) L^{n)^ 

(f>e^ 
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Obviously, L)(grad) = i?i(n) and thus, this dehnition coincides with the previous dehnition 
given in Example 11.21 Similarly, we dehne div as the closure of 

n 

i=l 

Moreover, we dehne the gradient on vector-helds, denoted by the same symbol grad, as the 
closure of 


and likewise the divergence of matrix-valued functions, again denoted by div, as the closure 
of 


i4’ik)i,ke{l,...,n} dk4>ik j 

\k=l / 

In the special case n = 3 we define the operator curl as the closure of 

Cooi^f C L^{Qf L^{Vlf 

f 4>l\ ( d2(t)3 - d'i(t)2 \ 

</>2 53(/)i - di ( j )^ . 

\4>3/ \di(p2- d24>l / 

By integration by parts we derive the following relations 

grad C —(div)* =: grad, 
div C —(grad)* =: div, 
curl C (curl)* =: curl. 

It is easy to see that for the scalar-valued gradient, we have Il(grad) = and thus, the 

dehnition of grad coincides with the previous one in Example 11.21 We recall that u G L)(grad) 
satishes n|aQ = 0, n G Zl(div) satishes n|ao • n = 0 and w G Zl(curl) satishes w\qq x n = 0 
for domains n with smooth boundary, where n denotes the unit outward normal vector held 
on on. As the operators grad, div and curl can be dehned for arbitrary open sets fl, we will 
use them to formulate the respective generalized boundary conditions, which do not need to 
make use of classical boundary traces. 


Beyond the aforementioned spatial differential operators, we need a suitable realization of the 
temporal derivative, see also |161 Example 2.3] or |ini Section 2] for a brief discussion. 
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Definition 2.2. Let ^) > 0 and H he a. Hilbert space. We consider the weighted L^-space 
H) := < f : M ^ H \ f measurable, f dt < oo 

[ Jr 

of iL-valued functions and equip it with the following inner product 

{f\9)gfi ■= j dt. 

R 

We define the derivative on as the closure of 

If the choice of g is clear from the context, we will omit the index g and just write do. 

Remark 2.3. (a) The subscript 0 in the operator do should remind of the “zero’th” coordinate, 
which in Physics’ literature often plays the role of the time-derivative. We shall adopt this 
custom here. 

(b) By definition the operator e~^ : Hgfl{^-,H) —)• {t is unitary. 

Consequently Cg := Te~^ : Hgfi{E.]H) —)• is unitary, where we denote by T : 

L‘^{W;H) —?• the unitary Fourier-transform given by 

(x) -=^1 e-“V(y) dy (X G M, </. G Coo(K; H)). 

R 

One can show that 

do = C*g{im + g)Cg, (2.1) 

where m : D{m) C L^(M;iL) —)• L^(M;Lf),/ i-)- (t i-)- tf{t)) with maximal domain, yielding 
that m is self-adjoint, and, by proving that do has is continuously invertible. Moreover, 

from (|2.1I) we read off that cr{do), the spectrum of do, coincides with i [M] -|- g. 

Let go > 0 and M : Be L{H) be analytic such that 

M{z)x\x)H > 0 (2.2) 

for every z e Be , x e H. Defining 

for g > go we can generate operator-valued multipliers and via M{dQ^) := £*M ^ i^+g ) 

corresponding operator-valued functions of We shall loosely refer to such operators 

as material law operators. 

We quote the following well-posedness result. 
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Theorem 2.4 f [161 Theorem 6.2.5], |151 Solution Theory]). Let H be a Hilbert space and 
A : D{A) C H ^ H a skew-self adjoint linear operator. Moreover, let > 0 and M : 

analytic such that there exists c > 0 with 

i)\e{z~^M{z)x\x)H > c|x||^ (2-3) 

for every z G Be ^ 2 ^ 1 ’ ^ ^ ^ H. Then for each Q > Qq the operator (j!)QM{df^) + 

is boundedly invertible on HgfliM.] H) and causal. 

Remark 2.5. (a) In many applications, M has the simple form M{z) = Mq + zMi for some 
Mq,Mi G L{H). Condition (12.31) can then for example be achieved for qq large enough, if Mq 
is selfadjoint and there exist cq, ci >0 such that 

{Mox\x)h > co\x\\ and y{t{Miy\y)H > ci\y\]j 

for every x belonging to Mq[H\, the range of Mq, and every y from N{Mq), the nullspace of 
Mq. The corresponding material law operator is then 

M{d^^) = Mq + 

(b) For the interested reader, we shall also mention possible generalization of Theorem 12.41 to 
non-autonomous equations mm) or non-linear equations (IMIEZ]) 


2.1 The Guyer-KrumhansI model of heat conduction. 

The Guyer-Krumhansl model of heat conduction (see e.g. 0015]) consists of two equations: 
A first equation, the balance equation, relates the heat 9 to the heat flux q in the way that 

gedoO + div q = h, 

where g and c are certain material parameters and h is a given source term. The difference 
to the classical equations of thermodynamics is the following alternative to Fourier’s law: 

(1 -I- tqOq) q = -K grad 9 giAq g ,2 grad div g, 

where tq, k, fii, g 2 are real numbers modeling the material’s properties. We mention that the 
choices tq = /ii = /r 2 = 0 recover the standard Fourier law. In any case, k is assumed to be 
strictly positive. With this observation, we may reformulate the modified Fourier’s law as 

(roK~^9o + — IjL 2 K~^ grad div) q = — grad0. 

Here we emphasize that here A - as does grad div - acts on vectors with 3 components. 
Summarizing, we end up with a system of the form 

/ gedo div = 

ygrad (roK“^9o -|-graddiv) J \q J J 

and thus, the block operator matrix to be studied reads as follows: 

/ gedo div 

Y grad [TQK~^dQ k~^ — g,iK~^A — g, 2 K~^ grad div) 
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In the following we show that the latter block operator matrix is of the form of abstract 
evolntionary eqnations as mentioned in the introdnction. Rearranging terms and separating 
the spatial derivatives from the time-derivative, the operator becomes 

^ / gc 0 \ / 0 0 \ f ^ ^ 

0 ToK~^ J yO J \grad (—graddiv) J 

In order to obtain easily accessible bonndary conditions leading to the well-posedness of the 
latter eqnation, we give a different representation of the term (—/riAc“^A —graddiv) 
on the basis of onr theory of abstract grad-div-systems. Beforehand, however, we will show 
the following lemma of combinatoric natnre: 

Lemma 2.6. Let fi 2 G K, G ]R>o with fii > —fi 2 - Then there exists C = C* & 

^(([^3x3, (j^3x3) gij-iQiiy positive definite such that for all (j) G Coo we 

(k“^/xiA -|- n~^g ,2 grad div) cj) = div C grad 4>. 


Proof. Withont restriction assnme that k = 1. As an ansatz, take 

C = aosyrriQ -|- aiP -|- a2skew (2-4) 

for ao,ai,a 2 G ]0,oo[ to be determined later on, where 

symo,sym,P,skew: ^ 


are defined by 


with trace : 


:= -trace*trace, 
3 


sym a := - (a + a 


skew a := 


T 


.T 


a — a 


symQ := (1 — P) sym = 
C, (T I—)• Observe that 


{a G , 

{a G , 

sym (1 — P) 


trace* : C ^ 

^ 1-^ 2;l3x3, 


/ 100 \ 

where 13 x 3 denotes the identity matrix in and, hence. Pa = I 0 10 1. Observing 

Vool/ 


^Recall that here A is the Laplacian acting as 


A 0 0 
0 A 0 
0 0 A 


and that div on the right-hand side is the row-wise 


divergence and the grad on the right-hand side is the Jacobian of (j). 
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that trace grad 4> = div 4>, we get that 
div syiiiQ grad = div 


di(l)j + djcjii - -6ij div( 




and 


dj^ (^di(j)j + dj(t)i - ^ div cj) 

= - A(/) + - grad div (j) — - grad div (j) 

Zt A o 

1 2 

= - (A(^ — grad div i?i) + - grad div 0 

Z o 

1 2 

= — - curl curl (j) + - grad div 4> 

Z o 


divPgrad 0 = divP 

'1 


ie{i,2,3} 


= div 


-6ij div (j) 


hie{i,2,3} 




*€{1,2,3} 


= - grad div cj). 

O 


Similarly, we compute 


div skew grad 4> = div skew {dk4>j)j 2 3 } 

= ^div(4</>i -5j4)j-fce{i,2,3} 


\\'^dk {dkcpj - dj(t)k) 


\k=l 


jG{l,2,3} 


= - (A — grad div) </> = -- curl curl (f). 

Thus, choosing C as in (12.41) . we get that 

div C grad 4> = div (aosyiug + aiP + a 2 skew) grad 4> 

= ao div symQ grad (f) + ai div P grad </> + 02 div skew grad cj) 

= ao ( — curl curl 4> -\— grad div 0 j + ai - grad div cj) — a 2 - curl curl cj) 
y 2 3 y 3 2 

ao + «2 1 , , , 2ao + ai 

curl curl cp -\ ---grad div cp 


2 

ao + a2 
2 

ao + a2 


(— curl curl + grad div) + 


2ao + ai ao + a2 


grad div 


ao + 2ai — 3 q;2 , ,. 

A H-grad div . 
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Comparing coefficients, we get 


OLQ + Ot 2 


= 


otQ + 2 q;i — 3(12 
6 


= ^^ 2 ■ 


Now, strict positive definiteness of C is equivalent ten ai, 02 , 0-3 > 0. Introducing A := 
oq — 3a2 £ IR) we get 

ai = 3 C 2 - ^ j , 

ao = l(m + 5), 


1 




A 


If /^2 > 0 then A = —3/xi < 0 is a possible choice for the latter set of equations to obtain strict 
positive definiteness of C. If /U 2 < 0 then A = 3(/i2—Aii) leads to /Ui + ^ = ^ 


0 , ^^ 2 - I =^^ 2 - ^ ^ 

which also implies that C is strictly positive definite. 


> 0 as well as /ri — 2 = /Ui — I {fj .2 — = 


_ 5/Ji-3At2 

2 


> 0 , 

□ 


Lemma 12.61 together with the previous remark provides a way for formulating the Guyer- 
Krumhansl model for heat conduction in the framework of evolutionary equations. This, in 
turn, results in the following well-posedness theorem, where we have for sake of definiteness 
chosen a specific case of boundary conditions: 

Theorem 2.7. Let tq,q,k G K> 0 ) Aii )/^2 £ 1^ with H 2 > —/J-i, 12 C open and C as in 
Lemma \ 2.fk Let c G L{L?‘{LL)) he selfadjoint and strictly positive definite and A: D{A) C 
L2 (L?)3 0 J2 0 j2 (q)3x3 ^ J2 (j^)3 ^ j2 (q) 0 j2 (^)3x3 

( 0 (grad — div) \ 

( 4 l) » 

for all 

( ffiA f L^nf \ 

G D{A) := M ^2 U L^n) 

iVV’3/ VL2 (0)3X3/ 

Then the (elosure of the) operator 

/ roK“3 0 0\ / 0 0 \ / 0 (grad -div) 

T := (9o ( 0 0 I + I 0 0 0 I + I / div A 

Vo OO/ \ 0 OC-^J VV -grad / 

is continuously invertible in L^ (0)3 0 L^ (0) © L^ (0)3x^^ for sufficiently large g > 

0 . 

^Note that P, symp and skew are projections on pairwise orthogonal subspaces of 
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Proof. First of all we observe that A is skew-selfadjoint. Indeed, ( ) : Hi (0)3 C 

V-grady 

(0)^ — 7 > L^( 0 ) 0 L^( 0 ) 3^3 generates an abstract grad-div-system (use Proposition II. Ql with 
Co = —grad, D{Co) = Hi{Ptf‘ and Ci = div) I and the skew-selfadjointness of A follows from 
Theorem 11.61 Therefore, the invertibility of T is guaranteed by the strict positive definiteness 
of c and C by Theorem 12.41 □ 


Thus, for h E Hg^o (M; (0)) and u = {ui,U 2 ,us) we have 



j-i 

0 0 0 
0 0 C-3 



0 (grad — div 



This is 


(tok ^doUi + K 0 grad — div ua = 0, 

gcdoU2 + divrti = h, 

— gradui = 0. 

Substituting - as a formal calculation - the third equation into the hrst, we get 
(tok^^SqUi + 0 grad n 2 — div C grad tti = 0 , 

which, by Lemma 12.61 gives the original form of the system back 

[ToK~^doUi 0 K~^ui) = K~^fj.iAui 0 «;“^/r 2 grad div til — gradn 2 . 


Remark 2.8. 


1. Note that the role of temperature is played by the second block component ■02 whereas 
ijji is the heat flux in this model. 

2. The system reveals appropriate possible generalizations for the material law operator. 
For example 

/ 0 0 \ / 0 0 \ 

M {do^) = 0 CO 0 0 5o"M 0 0 0 

Vo 00/ \ooc-V 

where k, kq, cq, C are a strictly positive definite operators acting in a matching (fl)- 
setting rather than just numbers. 


3. In proper tensorial terms the spatial operator should actually be written as 

0 ( div trace* — div) 

^ trace orad \ / 0 0 \ 


(trace grad \ 
-grad y 


0 0 
0 0 


0 

f trace \ 

V -1 J 


grad 


div (trace* 

0 0 


0 0 


-1) 


^As a subtlety not to be missed, here we have taken grad : H (P) ^ acting on scalars as the 

dual div^ of —div: —>■ This is in contrast to the situation in acoustics, where only 

the vanishing of the normal component on the boundary (or a generalization thereof) is imposed for the 
domain of —div. The tensor divergence occurring here is correspondingly div : H ^ , 

the dual of —grad on vector fields. 
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4. If preferred we may specialize to the symmetric tensor case, compare jlll 1^ . and use 

0 ( grad - div tsym ) 


div 

-i*ymgrad 


0 0 
0 0 


or, according to item 3, 


f trace \ * ° , 

I _1 1 ''symgrad 


div isym (trace* — 1 ) 


0 0 
0 0 


instead of the spatial operator above. Here we denote by igym : -^sym(^)^^^ 
the canonical embedding of 

{(/b W,2,3} G I (i,j € {1,2,3})} . 

In consequence, —)• is the orthogonal projection onto . 


2.2 A class of dynamic boundary conditions. 


In this section, we shall analyze a class of dynamic boundary conditions related to the wave 
equation. For this let H C be open and bounded with Lipschitz boundary F := and 
unit outward normal field n. Further, we shall assume that F is a manifold, where it is possible 
to define the co-variant derivative gradp as an operator from L'^iT) to L‘^(X) •= L"^ (r)^- By 
7 : H'i(H) L 2 (r),(/) !->■ (/)|r we denote the trace mapping. Consider the operator equation 


/ 


(9oMo + Ml + H) 


V 



\ 

/ 


for Q sufficiently large and with 


/ 

/ 

\ 


^ g ^ 



hi 


V' 

{h2/ 

V 


€ (M; l2(H) © l2(H) © L^{n) © lI{T)) 

(2.5) 


/ 


( 0 ) 


A = 


\ 


( grad 

7 

gradp 7 


( grad 

7 

gradp 7 


/O 0 0\ 
0 0 0 
\0 0 0 / 


\ 


/ 


where Mq is selfadjoint and strictly positive definite on its range, while Mi is strictly positive 
definite on the null space of Mq. Note that by the general solution theory Theorem 12.41 
well-posedness of the system just discussed is not an issue here, if we can ensure that A is 
skew-selfadjoint. For this, we shall have a closer look at the operator A. Defining for this 
section 


Xo := L^Q), 

Xi := (D(gradp 7 ), gradp 7 • Ii 2 (p) + I • ) , 

we realize that (|2.5I) can be treated as an abstract grad-div system: 
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Proposition 2.9. The operator 
( grad \ 

C := { 7 y.D (gradr 7 ) ^ Xq ^ © L\T) © lUt) 

\ gradr 7/ 

generates an abstract grad-div system. 

Proof. At first, note that Xi ^ Xq with dense and continuous embedding, yielding in par¬ 
ticular that C is densely defined. Moreover, as 7 : Hi{Q) —)• L‘^{T) is continuous and gradp : 
Hi{T) C L^(r) — L‘^{T) is closed, the operator gradp 7 : A>(gradp 7 ) C —)■ L‘^{T) is 

closed as well. Thus, C generates an abstract grad-div system by Proposition 11.91 for the 
choices Cq = gradp 7 , Ci = 7 and C 2 = grad. □ 


We shall now focus on the equations satisfied by 


tions, which are encoded in the domain of A. 



especially the boundary condi- 


Lemma 2.10. Let gradp 7 : Xi —)• (P). Then (gradp 7 )^ = 'j' gradp with gradp: Hi{T) —>• 

L^(r) and 7 : Xi —)• HifT) where 7 ': HifT)' —>■ X[ denotes the dual operator of 


Proof. For (/> G L^(r), we compute with if G (gradp 7 ) = Xi: 

((gradp 4>) (V’) = {(t>\ gradp 7 V’)L?©) 

= (gradp (f) {'yfj) 

= (Ygradpc/)) (V’). □ 

Remark 2.11. For divp := — gradf, we have the inclusion — 7 'divp C (gradp 7 )*. In fact, this 
follows from the previous lemma and the inclusion — divp C gradp, which holds by Lemma 

01 

Lemma 2.12. Let v G Zl(div), 71 G L‘^{T), 72 £ -f^i(r)', where HifT) = H(gradp) equipped 
with the graph norm. Consider the operators grad: Xi —>■ L^ (11)^ , 71 ; Xi —?• L^(r), xi eG 7 x 1 
and 72 : Xi —^ Hi{T),xi i-G 7 x 1 . Then —divx = grad^x + 7*71 + 7^72 if and only ijo 
n • X + 7 i + 72 = 0 on (F). 


Proof. (^) Assume that — divx = grad^x + 7 i 7 i + 7272 - Then, for G Xi we have that 
jijj = = 72-0 G F^i(r) G Hi/ 2 (T) G L^(r) and we compute 

(n - X + 71 + 72 ) ( 7 V’) = (divx|V^)p 2 (Q) + (x| grad- 0 )^ 2(!^)3 + (7il7i'0)L2(r) + 72 ( 7 V’) 

= (div x|V')p 2 (Q) + (grad^ x + 7^71 + 7372 ) (' 0 ) 

= (divx|V^)p2(Q) - (divx|V')p2(t2) 

= 0 . 

■‘Note that a-priori n-v is only defined as a functional in H 1/2 (^)^ where J7i/2(r) = 7 [ili(P)]. Indeed, n • v 
is the element on such that for all ‘ip £ ili(P) the equation 

(n • v) {yip) = (div'u|V')i 2 (n) + {x| gradi/))^2(n)3 

is satisfied. 
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Thus, the assertion follows, if we show that 7 [^i] is dense in For this, we note 

that C'oo(fl) is dense in Hi{^) since O is bounded and has Lipschitz boundary. Moreover, 
7 : Hi{Q) —>• Hi/ 2 (X) is continuous and onto and hence, 7 [ 6 * 00 (f^)] is dense in 
thus, dense in L^(r). Since C'oo(f^) C i 7 (gradp 7 ) = Xi we derive the assertion. 

(<^) Assume now that n • n + r/i + r /2 = 0 on (F). Then we compute for V’ £ Xi 
0 = {n-v + r]i+ r?2|7V')L2(r) 

= (divn|V’)L2(f2) + (n|gradV')^2(f2)3 + {vihii>)L^{r) + V2{l2i^) 

= (divn|V’)L2(n) + (grad^ n + 7^7/1 + 72772) W 

which gives the assertion. □ 


With this lemma we can characterize the domain of C* with C := 
abstract boundary condition. 

Theorem 2.13. ITe have 


( grad 

7 

gradp 7 



= — div V 


with 


/ grad \ 

I 7 I in terms of an 
\ gradp 7/ 


/ / grad \*\ 

T> I I 7 I I = {(77,771,772) G T) (div) X L^(r) X L^(r) I 77 • n + 771 + gradp 772 = 0 077 (F)} . 

\\gradr 7 / / 


Proof. In Proposition 12.91 we had already seen that C = 


grad 

7 I generates an abstract 
gradp 7 ^ 

grad-div-system and thus, by Theorem 11.61 we get (77,771,772) G D(C*) if and only if grad^n + 
7*771 + (gradp 7)^772 G L?‘{Pt). Moreover, we clearly have 


grad 

0 

0 


C 


grad 

7 

,gradp 7 _ 


and hence. Corollary 11.81 implies that 


C* C ( - div 0 0) , 

since — div = grad . In consequence (also use Lemma l2.10p . (77,771,772) G D{C*) if and only if 
77 G 11 (div) and — div v = grad* v + 7*771 + 7' gradp 772 G L^(n), which in turn is equivalent to 
77 G il(div) and 77 • 77 + 771 + gradp 772 = 0 on T^(r) by Lemma 12.121 □ 
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In the rest of this section, we shall formally compnte the eqnation modeled by (12.5p . That is 
to say, we assnme that the system 


(^oMo + Mi+A) 



/ 

/ 

\ 


( 9 ^ 


1 

hi 


V' 

Kh 2 j 



given in (|2.5I) with the special block-diagonal sitnation 


Mo + 


( moo 



(0 0 0 ) 

/ill 0 0 

0 /i22 0 
0 0 /i33 


+ <9o 


-1 


/ noo 



(0 0 0 ) \ 

liii 0 0 \ 

0 V22 0 j ’ 

0 0 Ji33 / / 


for non-negative scalars mooj noo, fin, fi 22 , fJ- 33 , ^^ 22 , z ^33 being arranged in the way that Mo 

and Ml satisfy the conditions to warrant well-posedness, i.e., Mq = Mq strictly positive 
definite on its range and IHeMi strictly positive definite on the nnllspace of Mq, admits 


smooth solntions 


f P \ 

(;) 

KWJ J 


The resnlting system leads to 


domoop + noop + divv = f 
dofiiiv + liuv + gradp = g 

dofi22Vi + ’^22Vi +lP = hi 
dofJ-33V2 + T^33m + gradp yp = /12 

n • n-I-? 7 i =-gradp ry 2 , (2.6) 

where the last eqnality is the bonndary condition indnced by the domain of A. Eliminating 
Vi ) ^2 by nsing the third and fonrth line of (12.6p , we obtain by formal calcnlations for the 
bonndary condition 

n-v + {dofi 22 + 2 ^ 22 )“^ {hi - 7 p) = - gradp (5o//33 + ^^ 33 )”^ {h 2 - gradp 7p) • 
Mnltiplication by (9o/^33 + ^^ 33 ) gives 

idofi33 + 1 ^ 33 ) n-v- {dofi 33 + 1233 ) idofi 22 + ^^ 22 )“^ 7P - gradp gradp 7P 

= - gradp /i 2 - {dofi33 + 1233 ) {doP 22 + 1 ^ 22 )"^ hi. 

With fi 22 = 0 and 1222 = 1 this simplifies fnrther to 

{dofi 33 + 1233 ) n-v - {dofi 33 + 1233 ) 7 p - gradp gradp 7 P = - gradp /12 - (9o/^33 + 1 ^ 33 ) h. 

Eliminating v, by nsing the second eqnation of (|2.6p . yields 

- {^o^J -33 + J 233 ) n ■ {dofiii + i 2 iiy^ {g - gradp) -k (5oP33 + 2 ^ 33 ) IP + gradp gradp 7 p 

= gradp /i 2 -h {dofi 33 + 1233 ) hi 
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and hence 


{dofJ -33 + 1 ^ 33 ) {doHii + i^ii) ^ n ■ gradp + { 80(133 + 1 ^ 33 ) 'yp + gradp gradp 7P 

= gradr ^2 + (9o^33 + 2 ^ 33 ) + { 80^33 + 1 ^ 33 ) {8oPii + vii)~^ n • g. 

Simplifying further by letting ^33 = pna, U 33 = vna with some non-zero a € M>o, we arrive 
at 


{8o(J‘ii + vii)'yp + n-gia.dp + a ^ gradp gradp 7P = (i9oA*ii + yii)hi + n- g + a ^gradr/ 12 , 

which is a heat type equation on the boundary of n (see also Remark 12. 111) . Noting that the 
first and second line of (|2.6I) gives 

8 omooP + nooP + div(ao/Uoo + t'oo)“^(5 - gradp) = /, 

which yields for vqo = 0 an abstract wave equation of the form 

8 lmooP + 8 onoQP - div poo gradp = 5o/ - div poo 9 - 

Thus, for suitable Mq and Mi, (12.51) models a wave equation on with a heat equation on 
90 as boundary condition. 

Remark 2.14. In |^, an equation was studied, where on a 1-codimensional submanifold in 0 
another partial differential equation on an interface is also considered. We shall outline here, 
how to include such problems in the abstract setting just discussed. So let O C open with 
Lipschitz boundary and assume that there exists O C 0 open with the property that 90 C 90 
and such that both 90 and 90 \ 90 are locally the graph of a Lipschitz continuous function 
and are manifolds allowing for the definition of a covariant derivative. As the operator C, 
generating the abstract grad-div system, we take the following operator 

/ grad \ 

C:= 7 

V gradr 7 / 

similar to the considerations above. The domain, however, changes a little, namely 7 being 
now the evaluation at T := 90 and gradp is given analogous to the above definition. 


2.3 A Leontovich type boundary condition. 


As for the acoustic equations just discussed, we formulate a similar problem also for Maxwell’s 
equations. We will show that Maxwell’s equations with impedance type boundary conditions 
can be formulated within the framework of evolutionary equations employing the notion of 
abstract grad-div-systems. For this let Q C R^ open. In this section, we consider the equation 





E (R; L^{Qf © © M^ace) • 


( 2 . 7 ) 


Here Mti-ace denotes a suitable Hilbert space, which will be used to formulate the boundary 
condition. In the forthcoming subsections we will discuss two possibilities for the choice of 

.^^trace • 
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2.3.1 A classical trace version 


Throughout, assume that is bounded and has Lipschitz-boundary T and denote by n the 
unit outward normal vector-held on T. In order to dehne A in (|2.7I) properly, we need the 
following trace operators (see m)- 

Definition 2.15. Let T^(r) denote the space of tangential vector-fields on T, i.e. 

Ll{T) :={/GL2(r)3|/.n = 0 }, 

which is a closed subspace of L^(r)^. We define the mappings 

tTt- : — >■ L^(r), LI I— —n x (n x yH) = yH — (n • 'yH)n, 

and 

y^ : Ll{T), H ^ x n, 

where y : —)• Hij 2 {^)^ denotes the classical trace (cp. Subsection 12.2p . We set 

Vt^ ■= TTr [ili(n)^] and := y,- [Lli(n)^], which are Hilbert spaces equipped with the norms 

Remark 2.16. 

(a) We have id 7 r: 14- ^ and id..^: ^ with continuous and dense embeddings. 

Consequently, id*: T^(r) as well as id*: T^(r) with continuous and dense 

embeddings. 

(b) Integration by parts gives 

(7r^H'|y^.E;)i2(r) = (curlLI|L;)i2(Q)3 - (il| curlL;)i2(Q)3, 
for all H,E ^ which yields that (see |2] or the concise presentation in |35l Section 

4 ]) 

o vr^ : H\nf C II(curl, Q) V', 

where R(^]^2y is defined as in Section[Tl is continuous, where II(curl, H) denotes the domain 
of curl equipped with its graph norm. Hence, has a unique continuous extension to an 
operator 

TTr : il(curl, H) Vj^. 

The same argument works for y,-, we, thus, get that 

'Jr : iI(curl,H) 


is continuous. 

Lemma 2.17. The operator 

C: vr-^ [id* [Ll (T)]] C L\nf ^ L^inf © L2(r), H ^ ~^ 

is closed. 
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Proof. First of all, note that vr^ ^ [id* [-L^(r)]] C F)(curl). Next, let {(j)n)n hr D{C) with 
4>n ^ 4> as n —)• oo and convergent to some ij: G L‘^(P) as well as {curl(j)n)n 

convergent to some rj G (O)^. By the closedness of curl, we infer that cj) G Ff(curl, n) 
and curl(/> = r]. Hence, (/>„ —5- </> in Ff(curl, H) as n ^ oo. As vr,- is continuous, we infer 
that TTT-cfn TTrf in as n —>■ OO. Since A^(r) with continuous embedding, we get 

TTr4> = Ip £ which implies the assertion. □ 

Dehne for this section 


Ao := 

X, := [d{C),V\-\^ + \C-\^). 


We may now introduce the operator A: 


/ /0 0 0 \ 
0 0 0 
\0 0 0 / 


A = 



\ 



\ 



/O 0 0\ 
0 0 0 
\0 0 0 / 
/O 0 0\ 
0 0 0 
\0 0 0 / 


\ 


/ / 


where curl := curl|xi : Ai ^ (Q)^ and := vr^-lxi : Ai ^ Lf(r). j ■ X^ C 

L^(Q)^ —)• 0 L^(r) generates an abstract grad-div system: the only thing left to show 

is that Xi is dense in L^(Q)^. This, however, is trivial as Coo(^) ^ X^. 


Similarly to the previous section, we shall have a look at the condition for being contained in 
—curl 


the domain of 


. Again, we need a prerequisite: 


Lemma 2.18. Let E G Il(curl), rj G L^(r). Then curlFl = curl E — 7 r*r/ if and only if 
'jrE + r] = 0 on L‘^ (T) . 


Proof. We observe that for 'I' G C Xi the equation 

{-/rE + rj) (vr^T) = (curlF;|T)^ 2(^)3 - (^FlIcurlT^^^^^^^ 0 (hl^r^)L 0 r) ( 2 - 8 ) 

= (curl.F|T)^2(^)3 - (curl^F; - (T) 

--O 

holds true. Thus, if curlFi = curl E — Trf rj, we get that {'jrE 0 rj) (vr,-'!') = 0 for each T G 
Thus, 'yT-E+rj = 0 on L^(r), due to the density of = 14- in A^(r), see [21 p 

850]. On the other hand, if Jt-E+t] = 0, equation (|2.8I) immediately gives curl FI = curl E—^P-rj 
by the density of in □ 
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Theorem 2.19. We have ( C ( — curl O) and 


D\\ j = {(T;,??) ET>(curl) xL2(r) | 7 ,£; + ?7 = 0 onL2(r)}. 


Proof. Note that with curl = curl* we have 


—curl\ /—curA _ ^ 

0 )-[ -iK 


Hence, by Corollary 11.81 we get 


C* C ( — curl 0) . 

Therefore, by Theorem 11.61 we obtain {E,ri) E D[C*) if and only if ill E ii(curl) and 

— curl if = —curl E + E 

which, in turn, by Lemma 12.181 is equivalent to 'jrE + r/ = 0 on (T) and E E Zl(curl). 


□ 


The latter theorem tells us that the containment in the domain of H is a boundary equation. 
So, in order to have a general class of Maxwell-type equations at hand, which include dif¬ 
ferential equations on the boundary, any suitable material law operator M{dQ ) allowing for 
well-posedness of (|2.7p is admitted. 

We shall elaborate the following particular choice for the material law operator: 


M ( 8 „-‘) 


/ h 

( 00 ) 

\ 



E L (k; L^nf ® A {nf © Ll (F))) 

\\v 

VO k(5o )) 

/ 


where k : Be i® analytic and satisfies condition (12.3p . For this choice 

(also use the result from Theorem 12.191) (|2.71) formally reads as 


dofiH + curl E = f, 
dosE — curlLf = g, 
doK (5(7^) r] + TwH = hi, 
•^rE + -q = f). 


Elimination of g yields 

— doK ^rE+ TTfH = hi- (2.9) 

This is an impedance type boundary condition. To see this we recall Jt-E = 'yE x n and 
tTt-H = — n X (n X yiL), use that (nx) (—nx) (nx) = (rex) and multiply (|2.9p by —rex. For 
ease of formulation, we shall suppress denoting the trace operator 7 . We obtain 

re X (Sore (Sf)"^) nxE)—nxH = —re x hi. 
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or, equivalently, 


—n X (c1ok(9q ^)n x {n x {n x -E))) = n x H — n x hi, 


which gives, 

TTrE = (^{—nx) doK {nx)) ^ {n x hi — H x n). (2.10) 

Now, in the literature, we find several choices for the impedance operator 

■= {{-nx) doK inx))~\ 

For discussing some of these choices, we let e', fj,', a' , t' G M>o be parameters with e' + n' > 0 
and + r' > 0. Mohsen |13j used the following form 

^ (^o"^) = {h' + r'd^^Y {s' + a'dYr^ > 


leading to a boundary condition, which Mohsen attributed to Rytov and Leontovich. Senior, 
|26) . and De Santis e.a., [3], discussed the case t' = 0, resulting in 

Z( 8 „- 

In the eddy current approximation {e' = 0) and also assuming r' = 0, a case discussed by 
Hiptmair, Lopez-Fernandez &: Paganini, [9], we arrive at the following form 






This leads to fractional evolution on the boundary. 

Remark 2.20. In |T] the author considers Maxwell’s equations in a two-dimensional setting 
and imposes a family of boundary conditions on a surface S contained in an exterior domain, 
which can be written in the form 


do YrE - TT^H) 





TTtH, 


( 2 . 11 ) 


where k denotes the curvature of S and 7 is an arbitrary function on S. It is shown that each 
of those boundary conditions yield a well-posed system and the asymptotic behavior of the 
corresponding solutions are studied. We remark here that the latter boundary condition is 
covered by (|2.9I) . Indeed, choosing 

(^0 + ^ + 7 ) 

and hi = 0 we obtain 

0 = -9oK (9(7^) JrE + TTrH 

= — ^1 — - ^(9o + 4 + "1^ ^ 

which gives 

0 = - ^«9o - ^ + 7^ IrE + (^^0 + ^ + 7^ 

which is just a reformulation of ( 12 . 111 ) . 
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2.3.2 An abstract trace version 


In this last section, we shall provide a possible way of by-passing boundary regularity require¬ 
ments giving the corresponding Leontovich type boundary conditions its most general form. 
For doing so, we need the following definitions, which are adapted from the abstract boundary 
data spaces, see e.g. [28l HSl [H] 

Definition 2.21. Recall that curl C curl, and consequently i/(curl, 11) is a closed subspace 
of Ff (curl, 12). We define the boundary data space o/curl by 

BV{cvlt\) := Fl(curl, 12)-*“, 

where the orthogonal complement is taken in //(curl, 12). Moreover, we denote by tcurl • 
BV{c\ir\) —)■ //(curl, 12) the canonical embedding, i.e. icuri'h = d*. Then, i*m.Fcuri = idBX)(curi) 
and 


Wi^uri : //(curl, 12) //(curl, 12) 

is the orthogonal projector on B'D{cxa:\). An easy computation shows that 

BV{c\\t\) = N{1 + curl curl), 

which in particular yields that curl[23/2(curl)] C 2322(curl). We set 

curl: BV{cxa:\) —^ .822(curl) 

^curl curl icurl 


Consequently, we have curlcurl= — idgx>(curi)- 


The main idea is now to replace the space T^(r) in the previous section by the space B'D{cui\) 
and the trace operator vr,- by the operator ^*^, 1 - Indeed, setting Xq := L^(12)^ and Xi := 

//(curl, 12) we obtain that ( : //(curl, 12) C L^(12)^ —)• L^(12)^ 0 .S/2(curl) generates 

V ^curl / 

an abstract grad-div system, see also Proposition 11.91 Consequently, (|2.7p with 


/ /0 0 0 \ 
0 0 0 
\0 0 0 / 


A := 


— curl 

''curl 


V 



\ 


/ 


V 


/o 0 o\ 
0 0 0 
\0 0 0 / 
/O 0 0\ 
0 0 0 
\0 0 0 / 


/O 0 0\ 
0 0 0 
\0 0 0 / 
/O 0 0\ 
0 0 0 
\0 0 0 / 




/ / 


is well-posed in //^^o L^(12)^ 0 L^(12)^ 0 BV(curl)'). In 

/_curl \ 0 

of I ^ ) ) we need to understand the operator 

V ''curl / 


order to characterize the domain 
: 23P(curl) —)■ //(curl, 12)'. 
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Lemma 2.22. TLe have = (1 + curP curl)tcuri- 


Proof. For ^ G H{curl,Q),^ G ;SP(curl) we compute 

= (icurl^|'I')H(curl,n) 

= (icurl^l^)L2(Q)3 + (curUcurl^l CUrl ^')i2(n)3 
= ((1 + curFcurl) tcurl^) (^'), 


which gives the assertion. 


□ 


Remark 2.23. Note that in Lemma 12.221 1 = R(^L 2 ^Q'^zy, the latter operator being given i 
Section [T] 

Thus, we end up with the following characterization result. 

Theorem 2.24. IFe have ( C ( — curl O) and 

V ^curl J ’ 

D ( ( 1 1 = \{E,r]) G L>(curl) x .SP(curl) 


m 


''curl 


^curl'^- curl 71 = 0 } . 


Proof. Since ^ ^ C ^ ^ , we obtain ^ ^ ^ — ( “ curl 0) by Corollary 11.81 

/ / - curl \ *\ 

Hence, by Theorem 11.61 we have {E,r]) G H ( I ^ only if Fi G Il(curl) and 

VV Curl J J 

— curl E = — curF E + rj 

= — curF Fi + (1 + curF curl)icurl^) 

by Lemma [2.221 The latter gives 

(curF — curl) Fi = (1 + curF curl) icurl^ 

= (— curl curl + curF curl)tcuri^ 

= (cUrF — curl) curl icurl^?; 

where we have used i3P(curl) = iV(l + curl curl). The latter is equivalent to 

CurlCurl curl ipm-iTy) 0. 

Indeed, we have that an element T G F)(curl) satisfies (curF — curl)'!' = 0 if and only if for 
each 4* G FI (curl) we have 

0 = ((curF — curl) T) ($) 

= (T| curl4>)i2(Q)3 - (curlT|4')i2(!^), 

which in turn is equivalent to '!' G Zl(curr) = F)(curl). Applying this to Fi — curlicurl^ £ 
F)(curl), we obtain the equivalence of 

(curF — curl) E = (curF — curl) curl tcurl?/ 
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and 


^curl^curl CUrl/-curl^) 0* 


( 2 . 12 ) 


Using the definition of curl and the injectivity of icurh we deduce the equivalence of (I2.12p to 

curl r] = 0. □ 

In the remaining part of this subsection, we like to point out that the boundary condition 
formulated in Theorem 12.241 is indeed a generalization of the classical boundary condition in 

Theorem 12.191 For doing so, let {E,r]) E T)(curl) x 5P(curl) with curl 77 = 0. Note 

• • • 

that the latter is equivalent to curl + 77 = 0, since curlcurl= —id^x'(curl)- Thus, we 

have to show that curl t*^j.iE can be identified with 'Jt-E in a certain sense. We compute for 
E,H e C°°{Ti) 

• • 

(curl ^x>{cur\) (^curl CUrl ^curl-^l^curhcurl'^)i^(curl,Q) 

= (curlicurl curl CUrUcurhcurl^)L2(n)3 + 

T (^curl curl 

~ (^curhcurl-®! CUrl (curl icurUcurl-® kcurFcurl-^)l/^ (H)® 

= -{E\ curli?)i2(Q)3 + (curlF;|i?)i2(s^)3 
= [{nx E)*H 


in case of a smooth boundary T := dVt. Thus, indeed nx E can be identified with curl 

or, in other words, curl: i3P(curl) —)■ .8P(curl) is a suitable generalization of the operator 

nx : T^(r) L^(r). Indeed, both operators are unitary and their adjoints are the negative 

• • 

operators, i.e. (nx)* = —nx as well as (curl)* = — curl . 
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